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Abstract: A non-singlet QCD analysis of the structure function X-F3 up to NNLO is 
performed based on the Bernstein polynomials approach. We use recently calculated NNLO 
anomalous dimension coefficients for the moments of the xF^ structure function in vN 
scattering. In the fitting procedure, Bernstein polynomial method is used to construct 
experimental moments from the xF$ data of the CCFR collaboration in the region of x 
which is inaccessible experimentally. We also consider Bernstein averages to obtain some 
unknown parameters which exist in the valence quark densities in a wide range of x and 
Q 2 . The results of valence quark distributions up to NNLO are in good agreement with the 
available theoretical models. In the analysis we determined the QCD-scale A.q^ d jy 4 = 
211 MeV (LO), 259 MeV (NLO) and 230 MeV (NNLO), corresponding to a s (M|) = 0.1291 
LO, a s (M|) = 0.1150 NLO and a s (M§) = 0.1142 NNLO. We compare our results for the 
QCD scale and the a s (M§) with those obtained from deep inelastic scattering processes. 
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1. Introduction 

The global parton analysis of deep inelastic scattering (DIS) and the related hard scattering 
data are generally performed at next-to leading order (NLO). Presently the next-to leading 
order is the standard approximation for most of the important processes in QCD. Analyzing 
DIS at next-to-next-to-leading order (NNLO) is important as we may be able to investigate 
the hierarchy LO — > NLO — > NNLO for the processes using the most precise available data. 

The corresponding one- and two-loop splitting functions have been known for a long 
time [1-11]. The NNLO corrections should be included in order to arrive at quantitatively 
reliable predictions for hard processes occurring at present and future high-energy colliders. 
These corrections are so far known only for the structure functions in the deep-inelastic 
scattering [12-15], for the Drell-Yan lepton-pair and gauge-boson production in proton- 
(anti-)proton collisions [16-19], and the related cross sections for Higgs production in the 
heavy-top-quark approximation [17,20-22]. 



Recently much effort has been invested in computing NNLO QCD corrections to a 
wide variety of partonic processes and therefore it is needed to generate parton distribu- 
tions also at NNLO, so that the theory can be applied in a consistent manner. Analysis on 
the NNLO cross sections for jet production is under way and it is expected to yield results 
in the near future, see Ref. [23] and references therein. For the corresponding three-loop 
splitting functions, on the other hand, only partial results have been obtained up to now, 
most notably on the lowest six/seven (even or odd) integer- N Mellin moments [24-26]. 
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S. Moch et al. [^7]] computed the higher order contributions up to three- loop splitting func- 
tions governing the evolution of unpolarized non-singlet quark densities in the perturbative 
QCD. 

During the recent years the interest to use CCFR data [28| for xF% structure function 
in the higher orders, based on the orthogonal polynomial expansion method has increased 
[29-34]. 

In this paper we determine the flavor non-singlet parton distribution functions, xu v (x, Q 2 ) 
and xd v (x, Q 2 ), using the Bernstein polynomial approach up to the NNLO level. This cal- 
culation is possible now, as the non-singlet anomalous dimension coefficients in iV-Moment 
space in three loops has already been introduced [^7], |3f| . 

The plan of the paper is to give an introduction to the CCFR data in Section 2. In 
Section 3 we present a brief review of QCD formalism of the non-singlet structure function 
in three loops. Parametrization of parton densities are written down in Section 4. Section 
5 contains a description of the Bernstein polynomial averages to be employed in the fits. 
Non-singlet quark distributions in the x-space are illustrated in Section 6. Section 7 con- 
tains a discussion and conclusions. 



2. CCFR experimental data 

The measurements of the CCFR collaboration provide a precise determination of the non- 
singlet deep inelastic scattering structure functions of neutrinos and anti-neutrinos on nu- 
cleons, xF^(x, Q 2 ). Data for xF% in neutrino-nucleon scattering is available from the CCFR 
collaboration |28|| . The data was obtained from the scattering of neutrinos off iron nuclei 
and the measurements span the ranges 1.26 < Q 2 < 199.5 GeV 2 and 0.015 < x < 0.75. 

The Q 2 -dependence of moments of structure functions can be predicted in perturbative 
QCD, and fits to data can be used to infer a s (M|). A difficulty is that there are upper and 
lower limits on the experimentally accessible range of x at low and high Q 2 , respectively. 
This is shown in Fig.l in which we plot the CCFR data for different values of Q 2 . We 
can see that at the lower range of Q 2 we are limited to low-x data, and at high Q 2 we are 
limited to the high x range. 

In order to reliably evaluate a moment at a particular Q 2 , we require data for the whole 
range of x. In fact for a given value of Q 2 , only a limited number of experimental points, 
covering a partial range of values of x are available. The method devised to deal with this 
situation is to take averages of the structure function weighted by suitable polynomials. 

Before reconstructing of the structure function from moments, we need to know how 
the structure function behaves in the missing data region. As we will see in next sections, 
in the fitting procedure we need to fit the xF% data of the CCFR collaboration. In this 
regard we can choose the extrapolation method. In this method, we fit xF^(x, Q 2 ) for each 
fixed value of Q 2 separately to the phenomenologically convenient expression 

xF^ hen \x) = Ax B {l-x) c . (2.1) 
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Figure 1: CCFR xF^ experimental data as a function of x and for some different values of Q 2 



This form ensures zero values for xF% at x = 0, and x = 1. The parameters A, B and C 
are obtained by performing x 2 fitting of Eq. fl2,l| ) to data for They are Q 2 -dependent 
quantities, and errors on their values are obtained by performing the fitting with the data 
for XF3 shifted to the two extremes of the error bars. 

In Table 1 we have presented the numerical values of A,B and C at Q 2 = 20, 31.6, 
50.1, 79.4, 125.9 GeV 2 . We have only included data for Q 2 > 20 GeV 2 , this has the merit 
of simplifying the analysis by avoiding evolution through flavor thresholds. 
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Q 2 (GeV 2 ) 


A 


B 


C 


20 


4.7425 ± 0.7585 


0.6356 ± 0.0638 


3.3756 ± 0.0892 


31.6 


5.4735 ± 0.8722 


0.6936 ± 0.0456 


3.6595 ± 0.0537 


50.1 


5.6795 ± 0.4841 


0.6977 ± 0.0843 


3.8390 ± 0.0459 


79.4 


4.5076 ± 0.5033 


0.5666 ± 0.0808 


3.7574 ± 0.0486 


125.9 


7.0775 ± 0.8641 


0.8186 ± 0.0475 


4.2456 ± 0.0458 



Table 1: Numerical values of fitting A,B,C parameters in Eq. ( |2.1|) . 



3. QCD formalism 

To carry out the analogous analysis at NNLO we need both the relevant splitting functions 
as well as the coefficient functions. Now not only the deep inelastic coefficient functions are 
known at NNLO, but also the anomalous dimensions in iV-Moment space are available at 
this order [27, 35 1. In this section we want to introduce the non-singlet structure function 
in Mellin moment space up to three loops order. 

The structure function xF%, associated with the parity- violating weak interaction, rep- 
resents the momentum density of valence quarks. So in the LO approximation we can 
write, 

xF% N = xu v (x) + xd v (x) + 2 xs(x) — 2 xc(x) , (3-1) 
xF$ N = xu v (x) + xd v (x) — 2 xs(x) + 2 xc(x) , 

where u v = u — u and d v = d — d are the proton valence densities. The asymmetry of the 
s — c doublet results in xF^ N xF^ N . The method of extracting CCFR experimental 
data, extracts the average of the neutrino and anti-neutrino distributions, so that 

xF 3 (x) = — — ^ Xi?3 = xu v (x) + xd v (x) , (3.2) 

here it is obvious that the xF^(x) is related to the combination of valence quark densities. 
Let us now define the Mellin moments for the vN structure function xF^(x, Q 2 ): 

M V 3 N (N, Q 2 ) = f x^F^x, Q 2 )dx . (3.3) 
Jo 

The theoretical expression for these moments obey the following renormalization group 
equation [29| 

^ + (3{As) jA s +7n S (A s )\M(N,Q 2 /» 2 ,A s (^)) = . (3.4) 

The symbol A s denotes the strong coupling constant normalized to A s = ot, s /{Air) and is 
governed by the QCD /^-function as 

^=(3(A s ) = -2Y,AAl +2 . (3.5) 
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c 



Eq. ( |3.5| ) is solved in the MS-scheme applying the matching of flavor thresholds at Q 2 = m 
and Q 2 = m\ with m c = 1.5 GeV and nib = 4.5 GeV as described in [36, 37]. MS- 
scheme convention introduced in 35 ] is extended in this way. In order to be able to make a 
comparison with the other measurements of Aqcd we adopt this prescription. The solution 



of Eq.flUD in the NNLO is given by 



1 Aln(lnQ 2 /Af^) 



^(l nQ 27A^3 $ ln '( ln Q 2 / A lls) ~ # H^QV^Ts) + ft/% ~ P\] ■ 



MS' 



(3.6) 



Notice that in the above the numerical expressions for Pq, Pi and P2 are 

p = 11 - 0.6667/ , 
Pi = 102 - 12.6667/ , 

P 2 = 1428.50 - 279.611/ + 6.01852/ 2 , (3.7) 
where / denotes the number of active flavors. 



In Mellin-iV space the evolution equation is solved |32j. The nonsinglet structure 
function M. (N, Q 2 ) is given by 

M{N,Q 2 ) = [1 + C^(N)A S + C^{N)A 2 ] f NS (N, Q 2 ) , (3.8) 

where the f NS (N,Q 2 ) is the Mellin transform of the non-singlet quark combinations and 
the C( k \N) are the corresponding Wilson coefficients 3£j. For the remainder of this paper 



we simplify our notation by dropping the sub- and superscript l vN' and '(3)' in Eqs. J3.4j , 
The solution of the non-singlet evolution for the parton densities to 3- loop order reads 



/*W) = f NS (N,Ql) ( ^g) 7 ° m {l + iMQ 2 ) - M Q 2 o)] (f-f)| 



+ [A S (Q 2 )-A S (Q 2 )] 2 ^{^ 



Pl(j 1 _jl^_ 
8P 2 \Pi Po 



+\[Al{Q") - A 2 M)} (jl. ~ |ti + £=^tf S ) \ , (3.9) 
where f NS is the valence quark compositions as 



f NS = { u -u) + (d-d) . (3.10) 
By considering symmetry between sea quark distributions we can write 

f NS (N, Ql) = u v (N, Q 2 ) + d v (N, Ql) . (3.11) 
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In the next section we will introduce the functional form of the valence quark dis- 
tributions and we will parameterize these distributions at the scale of Qq. As we see in 
Mellin-iV space the non-singlet parts of structure function in the NNLO approximation, 
i.e. A4(N,Q 2 ), can be obtained from the corresponding Wilson coefficients C^ k \N) and 
the non-singlet quark densities. 

By using the anomalous dimensions in one, two and three loops from [27] and inserting 
them in Eq. ( |3,9| ) and using Eq. ( |3.8| ), the moments of non-singlet structure function as a 
function of N and Q 2 are available. The results of [27] for 7 JVS+ (n) are depicted in Fig. || 
for four active flavors and typical values a s = 0.15,0.2 for the strong coupling constant. 




4. Parametrization of the parton densities 

In this section we discuss how we can determine the valence quark densities at the input 
scale of Qq = 1 GeV 2 . First of all, we should notice that the sum of the u v and d v 
distribution functions can be obtained from the CCFR data and not the two distributions 
separately. To separate the xu v and xd v contributions to xF% we need to use the relation 
between two distribution functions. 

To start the parameterizations of the above mentioned parton distributions at the 
input scale of Qq, we choose the following parametrization for the u- valence quark density 

xu v (x,Ql) = N u x a (l-x) b (l + c^ + dx) . (4.1) 

In the above the x a term controls the low-x behavior parton densities, and the (1 — x) b 
terms the large x values. The remaining polynomial factor accounts for additional medium- 
x values. To separate the xu v and xd v contributions to xF^ we assume the relation between 
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two distribution functions as 

d v 



M{1 - x) e , (4.2) 



this equation is same as the ratio which has reported in Ref. p0|] . The parameter e in this 
ratio is very important to control the behavior of F3 for large x value and the coefficient M 
is the normalization constant {= jf^)- Therefor the parametrization of the d- valence quark 
density is as follows 



xd v (x, Qq) = — -(1 - x) e xu v (x, Qq) . (4.3) 
Normalization constants N u and Na are fixed by 

f u v (x)dx = C u , (4.4) 
J 

-1 

d v (x)dx = C d , (4.5) 





so normalization constants are equal to 

N. = — (A 6) 

B(a, 1 + b) + aB(l/2 + a, 1 + b) + dB(l + a, 1 + b) ' v ' ' 

d B(a, 1 + b + e) + cB(l/2 + a, 1 + 6 + e) + dB(l + a, 1 + 6 + e) ' 1 ' j 



here C u = 2 and = 1 are respectively the number of u v and quarks and B(a,b) = 
r(a)r(b) j g ^ e g u j er g e ^ a function. 

1 (a+b) 

The above normalizations are very effective to control unknown parameters in Eqs. ( [4.1| , f4.3| ) 
via the fitting procedure. The five parameters with Aq£^ will be extracted by using the 
Bernstain polynomials approach. 



Using the valence quark distribution functions, the moments of u v (x, Qq) and d v (x, Qq) 
distributions can be easily calculated. The Mellin moments for the sum of the two valence 
quark distributions in the proton is as follows 

u v (N,Ql)+d v (N,Ql) = [ x N ~ 2 (xu v (x,Ql) + xd v (x,Ql)) dx 

Jo 
-1 

JV-2 



X 







xu v (x,Ql) ^1 + ^L(1- X )^ dx. (4.8) 



Now by inserting the above equation in the Eq. ( 3.1 1| ) , the function of / ns (N,Qq) is 



determined in terms of unknown parameters a, b, c, d, e. This function is needed to 
determine the moments of non-singlet structure function in the related order. 
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5. Reconstruction of the structure function from moments 



Although it is relatively easy to compute the iVth moment from the structure functions, 
the inverse process is not obvious. To do this, we adopt a mathematically rigorous but easy 



method [41| to invert the moments and retrieve the structure functions. The method is 
based on the fact that for a given value of Q 2 , only a limited number of experimental points, 
covering a partial range of values of x are available. The method devised to deal with this 
situation is to take averages of the structure function weighted by suitable polynomials. 
We define the Bernstein polynomials as follows, 



B 



nk (.•&) 



Tin + 2) 



T(k + l)r(n - k + 1 



-x k (l - x) 



n—k 



n> k. 



(5.1) 



These polynomials have a number of useful properties. These functions are normalized such 
that Jq B n ^(x)dx = 1 and they are also constructed such that they are zero at endpoints 
x = and x = 1. These polynomials are positive and have a single maximum located at 



%nk\%) 



x B nk {x) dx 



T(n + 2) T{k + 2) 



T(n + 3) r(fc + 1) 

and finally, they are concentrated around this point, with a spread of 



(5.2) 



(x - x nk ) B nk (x) dx 



T(n + 2)r(fc + 3) /T(n + 2) T(k + 2) 



V T(n + 4) T(A; + 1) \T(n + 3) T(k + 1) 
Therefore, for a given value of Q 2 , the Bernstein averages of F% which are defined by, 



(5.3) 



F 



nk 



(Q 2 )= f 

Jo 



dxB nk (x)F 3 (x,Q 2 



(5.4) 



represents an average of the function Fs(x, Q 2 ) in the region x nk — \^x nk <x<x nk + \lS.x n k- 
The key point is, the values of F% outside this interval have a small contribution to the 
above integral, as B nk {x) tends to zero very quickly. By a suitable choice of n, k we 
manage to adjust to the region where the average is peaked around values which we have 



experimental data |2fi 



The construction of an acceptable average, and the resulting suppression of the missing 
data region is demonstrated in Fig. |3[ In this figure the light grey region represents the 
interval x nk — \&-x n k<x<x n k + ^Ax n k and the dark grey areas represent the missing data 
regions. The small size of the dark grey region in the right hand plot demonstrates that 
this average has a negligible dependence on the missing data regions. Note that the right 
hand plot actually shows the integrand of the Bernstein average. The average itself will be 
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this function integrated over [0, 1]. 



By expanding the integrand of Eq. ( p~^ ) in powers of x, we can relate the averages 
directly to the moments, 



and using the definition of Mellin moments of any hadron structure function, Eq. ( |3.3| ), we 
have 

F„*(Q*) - £ jjj^^ljj ««* + + l.Q 2 ) • <M> 

We can only include a Bernstein average, F n k, if we have experimental points covering 
the whole range [x n k — ^Ax n k,x n k + \A.x n k] [ 42| ,[ 43] , [44 ] . This means that with the 



available experimental data we can only use the following 28 averages, including f!^ p \q 2 ) 



F 3 ( f p) (Q 2 ), F^ xp) (Q 2 ), ... . Using Eq. the 28 Bernstein averages F nk {Q' 2 ) can be 



written in terms of odd and even moments. For instance: 

F 21 (Q 2 ) = 6(M(2,Q 2 )-M(3,Q 2 )) , 

F U (Q 2 ) = 24 (0.5 M(2,Q 2 ) - M(3,Q 2 ) + 0.5 M(A,Q 2 )) , 

F&(Q 2 ) = 60 (0.5 M(3,Q 2 ) - M(4,Q 2 ) + 0.5 M(5,Q 2 )) , 

i (5-7) 

We can now compare the theoretical predictions with the experimental results for the 
Bernstein averages. Another restriction we assume here, is to ignore the effects of moments 
with high order n which do not strongly constrain the fits. To obtain these experimental 
averages from CCFR data p3], we fit xF^(x, Q 2 ) for each bin in Q 2 separately to the phe- 



nomenologically convenient expression given in Eq. (|2.1[). Using Eq. with the fitted 



values of A, B and C, one can then compute F^ p \q 2 ) using Eq. (|5.5|) , in terms of Gamma 



functions. Some sample experimental Bernstein averages are plotted in Fig. |I]in the higher 
approximations. The errors in the ^it (Q 2 ) correspond to allowing the CCFR data for 
X-F3 to vary within the experimental error bars, including the experimental systematic and 
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statistical errors [28]. 

The unknown parameters according to Eqs. ( fP| , |4.3| ) will be a,b,c,d,e and Aq^ d . Thus, 
there are 6 parameters for each order to be simultaneously fitted to the experimental 
FnkiQ 2 ) averages. Using the CERN subroutine MINUIT p5|j , we defined a global x 2 
for all the experimental data points and found an acceptable fit with minimum x 2 /dof = 
92.259/134 = 0.688 in the LO case, 77.452/134 = 0.578 in the NLO and 74.772/134 = 0.558 
in the NNLO case. The best fit is indicated by some sample curves in Fig. The fitting 
parameters with their uncertainties and the minimum x 2 values in each order are listed 
in Table 2. From Eqs. (4. 1,4~3|), we are now able to determine the xu v and xd v at the 




Figure 4: NLO and NNLO fits to Bernstein averages of xF^. 



scale of Qq in higher order corrections. In Fig. ^ we have plotted the NLO and NNLO 
approximation results of xu v and xd v with correlated errors at the input scale Qq = 1.0 
GeV 2 (solid line) compared to results obtained from NNLO analysis (left panels) and NLO 
analysis (right panels) by BBG ||| (dashed line), MRST (dashed-dotted line) J!?]] and 
A05(dashed-dotted-dotted line) @. 
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L(J 


ATT f ~\ 


ATfvTT (~\ 
1N1NLU 


u v 




1.952 


3.942 


5.134 




a 


U.57U ± 0.01b 


0.777 ± 0.020 


0.830 ± 0.019 




b 


3.300 ± 0.089 


3.548 ± 0.032 


3.724 ± 0.039 




c 


-0.380 ± 0.040 


0.410 ± 0.039 


0.040 ± 0.004 




d 


4.900 ± 0.136 


1.500 ± 0.027 


1.449 ±0.104 




N d 


1.308 


2.620 


3.348 




e 


1.699 ±0.135 


1.563 ±0.014 


1.460 ±0.020 




211 ±27 


259 ± 18 


230 ± 12 


X 2 /ndf 


92.259/134 = 0.688 


77.452/134 = 0.578 


74.772/134 = 0.558 



Table 2: Parameter values of the LO, NLO, and NNLO non-singlet QCD fit at Q% = 1 GeV 2 . 




-3 -2 -1 v -1 -2 -1 

10 10" 10 x 10" 10" 10 



Figure 5: The parton densities xu v and xd v at the input scale Qq = 1.0 GeV 2 (solid line) compared to 
results obtained from NNLO analysis (left panels) and NLO analysis (right panels) by BBG |^6| (dashed 
line), MRST (dashed-dotted line) @ and A05 (dashed-dotted-dotted line) @. 

All of the non-singlet parton distribution functions in moment space for any order 
are now available, so we can use the inverse Mellin technic to obtain the Q 2 evolution of 
valance quark distributions which will be done in the next section. 
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6. Valence quark densities in the x-space 



In the previous section we parameterized the non-singlet parton distribution functions at 
input scale of Qq = 1 GeV 2 in the LO, NLO and NNLO approximations by using Bernstein 
averages method. To obtain the non-singlet parton distribution functions in x-space and 
for Q 2 > Qq GeV 2 we need to use the non-singlet evolution equation for parton densties 
to 3- loop order in Eq.( |3.9D . To obtain the x-dependence of parton distributions from the 
./V— dependent exact analytical solutions in the Mellin-moment space, one has to perform 
a numerical integral in order to invert the Mellin-transformation 



1 f°° 'ik 
q v (x, Q 2 ) = - / dwIm[e^x- c - we ^M qv (N = c + we**, Q 2 )] , (6.1) 
71" Jo 

with q v = u v ,d v . In this equation the contour of the integration lies on the right of all 
singularities of M qv (N = c + we 1 *, Q 2 ) in the complex iV-plane. For all practical purposes 
one may choose c ~ l,<p = 135° and an upper limit of integration, for any Q 2 , of about 
5 + 10/ In x -1 j instead of infinity, which guarantees stable numerical results ||50| , |5lf| . In 
this way, we can obtain all valence distribution functions in fixed Q 2 and in x-space. In 
Fig. |6]we have presented the parton distribution xu v at some different values of Q 2 . These 
distributions were compared to LO , NLO and NNLO approximations with some theoretical 
models g7|, g|, ||. 

In Fig. |7] we have presented the same distributions for xd v . We should notice that in 
Figs. ||, H and ^ the minimum value of Q 2 from [47] and [52| is 1.25 GeV 2 and 1.3 GeV 2 
respectively. 

In Table 4 comparison of low order moments at Q 2 = 4 GeV 2 from our non-singlet NNLO 
QCD analysis with the NNLO analysis BBG06 |4|], MRST04 @, A02 @ and A06 §3] 
has been done. 



/ 


N 


NNLO 


BBG06 


MRST04 


A02 


A06 


u v 


2 


0.2934 ± 0.0036 


0.2986 ± 0.0029 


0.285 


0.304 


0.2947 




3 


0.0825 ± 0.0012 


0.0871 ± 0.0011 


0.082 


0.087 


0.0843 




4 


0.0311 ± 0.0004 


0.0333 ± 0.0005 


0.032 


0.033 


0.0319 


d v 


2 


0.1143 ± 0.0013 


0.1239 ± 0.0026 


0.115 


0.120 


0.1129 




3 


0.0262 ± 0.0003 


0.0315 ± 0.0008 


0.028 


0.028 


0.0275 




4 


0.0083 ± 0.0001 


0.0105 ± 0.0004 


0.009 


0.010 


0.0092 



Table 4: Comparison of low order moments at Q 2 = 4 GeV 2 from our non-singlet NNLO 
QCD analysis with the NNLO analysis BBG06 @, MRST04 @, A02 |4§] and A06 || 



7. Conclusion 



The QCD analysis is performed in LO, NLO and NNLO based on Bernstein polynomial 
approach. We determine the valence quark densities in a wide range of x and Q 2 . The 
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Figure 6: The parton distribution xu v at some different values of Q 2 . The solid line is our model, 
dashed line is the MRST model @, dashed-dotted line is the A05 model @ and dashed-dotted-dotted 
line is the CTEQ model pi. 



QCD scale Aq^^ is determined together with the parameters of the parton distributions. 



In Table 5 we have summarized our fit results comparing Aq^^ and a s (M§) for the LO, 
NLO and NNLO analysis. The LO value of Aq £ c = 4 is found to be smaller than the NLO 
value, while the NLO value comes out somewhat higher than the NNLO value. 





Aq^d, MeV 


a s (Ml) 


LO 
NLO 
NNLO 


211 ±27 
259 ± 18 
230 ± 12 


0.1291 ±0.0025 
0.1150 ±0.0011 
0.1142 ±0.0008 



Table 5: Aq^ and a s (M§) at LO, NLO and NNLO. 



We compare the results of the present analysis to results [34], [46|, p8| , [52-57] obtained 
in the literature at NLO and NNLO in Table 6, where most of the NLO values for a s (M|) 
presented are determined in combined singlet- and non-singlet analysis. The NLO values 
for a s {M%) are larger than those at NNLO in several analysis. The difference of both 
values, however, is not always the same. This is most likely due to the type of the analysis 
being performed (singlet and non-singlet, non-singlet only, etc.), in which also partly dif- 
ferent data sets are analyzed. Non-singlet QCD analysis were also performed for neutrino 
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Figure 7: The parton distribution xd v at some different values of Q 2 . The solid line is our model, 
dashed line is the MRST model @, dashed-dotted line is the A05 model @ and dashed-dotted-dotted 
line is the CTEQ model fl. 



data by using Jacobi polynomial method. In [33] the CCFR iron data on xF 3 (x,Q 2 ) |§ 
were analyzed in NLO and NNLO using fixed moments. Likewise a NNLO analysis was 
performed in |34|. In 33] rather large values for Aq^' MS : Aq^'^q = 371 ± 72MeV, 

Aqcd nnlo = ^16 i 51 MeV are obtained, which are larger than the values obtained in 
the analysis based on F% ,d (x, Q 2 ) data, still showing the pattern that the NNLO value 
is lower than the NLO value. In Ref. @ one finds AjJgJJK) = 281 ± 57 MeV, 
A Sc = D;NfLO = 2 55±55MeV. 

And finally in [^] with the QCD analysis of deep inelastic world data, the value of 
AqcI) ' MS * s reported as 

Aqc^S) = 265 ±27MeV (7.1) 
Aqct'Nf lo = 22 6 ± 25 MeV . (7.2) 

which seems close to results of the present analysis. We believe that the difference of the 
reported value above, not only depends on the type of analysis being performed (singlet 
and non-singlet, non-singlet only, etc.) but also on the kind of approach (iV-space, x-space, 
etc.) have been taken. 
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a s (ivi z ) 


Rpf 


NLO 






CTEQ6 


0.1165 ± 0.0065 


n 


A AT> GTH9 

MKb 1 03 


nil r* c i n nnon 

0.1165 ± 0.0020 


[fT/ili 




n 1 1 7i -i- n nm ^ 


U±gJ 


ZEUS 


1166 ± 0049 


IMMI 


HI 

111 


U.llOU ZE U.UU1 1 


IPhI 


\jl\jO 


119 


ipju 


BBG 


1148 ± 0019 


[Ell 




n 1 1 -u n nm i 

U.llOU HZ U.UU11 




NNLO 






MRST03 


0.1153 ±0.0020 


m 


A no 

A02 


0.1143 ± 0.0014 


IXcIl 


SYOlfen"! 


1166 ± 0013 




SYOl(^N) 


0.1153 ±0.0063 





GRS 


0.111 




A06 


0.1128 ±0.0015 




BBG 


11 34+0-0019 
U.110^_q 0021 


a 


Model 


0.1142 ±0.0008 





Table 6: Comparison of a s (M|) values from NLO, and NNLO QCD analysis. 

Another important characteristic of the deep inelastic neutrino-nucleon scattering is 



the Gross-Llewellyn Smith (GLS) sum rule [58] 

GLS(Q 2 ) = ■ 



xF* p+vp 



(x,Q 2 



-dx . 



(7.3) 



In the work of Ref. |]59| ], the following result of the measurement of the GLS sum at 
the scale | Q 2 |= 3 GeV 2 was reported: 

GLS(\ Q 2 |= 3 GeV 2 ) = 2.5 ± 0.018(stat.) ± 0.078(syst.) (7.4) 

In Ref. |6(| the GLS(| Qq |= 3 GeV 2 ) were analyzed based on the Jacobi polynomials 
expansion method. The value of GLS in the NLO approximation is reported as GLS (3 
GeV 2 )=2.446±0.081 @, which is in agreement with the results Eq. ( |7.4| ), We should 
notice that in order to obtain NLO expression for the GLS sum rule one should consider 
the NNLO approximation of the moments M(N, Q 2 ). Following Ref. we analyze the 
GLS sum rule with the corresponding perturbative QCD predictions for the first Mellin 
moments, and obtain 

GLS(\ Q 2 |= 3 GeV 2 ) = 2.40 ± 0.06 (7.5) 

We hope to report on the application of the methods employed in the present work 
to describe more complicated hadron structure functions, and on using the singlet case to 
extract parton densities in three loop in future works. 
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